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ABSTRACT 

Tidal radii of remote globular clusters {Roc ^ 35 kpc) are used to provide constraints 
of the mass profile of the Milky Way galaxy that are independent of kinematic data. 
The available data are consistent with the profile of an isothermal sphere with circular 
velocity Vc — 220 ± 40 km/s in the radial range 35 kpc ^ Rgc ^ 100 kpc, in good 
agreement with all recent estimates. The more robust constraint at large distances from 
the galactic center is provided by NGC 2419, yelding an enclosed mass of 1.3 



+2.9 
-1.0 



10^^ Mq at Rgc - 90 kpc. 
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1 INTRODUCTION 

One of the most straightforward examples of the great difh- 
culties associated with the measure of basic physical quanti- 
ties on the astrophysical scales is provided by the quest for 
the mass of the Milky Way. Even applying the most refined 
analysis and using the whole wealth of available data, real- 
istic uncertainties affecting single estimates of such a funda- 
mental parameter tipycally amount to ~ 100 — 300 % (see, 
e.g. .Wilkin son fc Evans 1999, for a state-of-the-art analysis). 
The whol e problem was r ecently reviewed and critically dis- 
cussed bv lZaritskvl il999l hereafter Z99). This author notes 
that the concept of total mass of the Milky Way is some- 
how ill-defined since we ignore the actual extent of the Dark 
Matter (DM) halo of the Galaxy. Hence, it is much safer to 
refer mass estimates to the enclosed mass within the galac- 
tocentric distance (Roc) sampled by the mass-tracer under 
consideration. This approach allows a sensible comparison 
between different estimates, since consistency requires that 
all estimates shall be in agreement with a unique mass profile 
[M{Rgc)] over the whole range of galactocentric distances 
that can be probed. Z99 uses the isothermal sphere as a ref- 
erence model to perform such a comparison and concludes 
that all the available estimates, covering different ranges in 
Rgc and using different tracers (from the HI rotation curve 
to the outermost Galactic satellites) , are consistent with the 
mass profile of an isothermal sphere with rotational velocity 
Vc ~ 180 km/s. Note that in this context the isothermal 
sphere is (obviously) not intended as a realistic model for 
the Galactic DM halo, but just as a suitable mass profile to 
provide comparison and cross-validation of the various es- 



E-mail: bellazzini@bo.astro.it 



timates. Hence, despite the large uncertainties affecting the 
single mass estimates, a remarkable overall consistency is 
apparent. If interpreted in the classical Newtonian frame- 
work, the mass of the Galaxy (Mc) is observed to grow 
approximately as Mq oc Rgc and the enclosed mass within 
Rgc ^ 300 kpc is of order Mg ~ 2 x 10^^ Mq (Z99). 

All the estimates considered by Z99 (as well as the more 
recent ones in the post- '99 literature) rely on the kinematics 
of the adopted tracers (rotational velocity curve, escape ve- 
locity of local stars, motion of satellites, either stars, glob- 
ular clusters or dwarf galaxies). Therefore, an estimate of 
Mg{Rgc) not based on kinematical data would provide 
a further important consistency check of our ideas on the 
mass profile of the Galaxy. Such kind of probe may be pro- 
vided by tidal radii of globular clusters (GC) ivon Hoernerl 
11959: iKina . 19621 Theory predicts that the Galactic tidal 
field fixes the cut-off in the density profile of GCs, the po- 
sition of the cut-off depending on the cluster distance and 
on the ratio between the mass of the Galaxy and the mass 
of the cluster (rric). Hence, having and Rgc from ob- 
servations, an estimate of Mg may be obtained, though af- 
fected by large uncerta inties. The approa ch has been at- 
tempted in the past by 'Wakamatsu' fl98l', hereafter W81) 
and by Innancn. Harris & Wcbbink (1983, hereafter IHW). 
W81 used globular clusters and dwarf spheroidal galaxies as 
tracers, and obtained Mg (78 kpc) = 'i.Otlir x lO" Mq. 
The result is in marginal disagreement with the conclu- 
sions by Z99 but is probably affected by the large uncer- 
tainties in the distance and tidal radii of dwarf spheroidals 
(for the same reasons these galaxies are not considered 
in the present analysis). IHW used only GCs and found 
Mg (44 kpc) = 8.9 ± 2.6 x lO" Mq, consistent, within the 
uncertainties, with the Z99 results. 
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There are cogent reasons to try to repeat the ex- 
periment twenty years later. The quahty of the ob- 
servational material is greatly improved, mainly thanks 
to the extensive compilation of surface brightness pro- 
files by iTraeer. King fc Diorgovskil lll995h . and theoret- 
ical advancements put t he pr o blem in a new light 
JOh. Lin fc AarsethI Il99l llOQ^ iMeziane fc ColinI Il996l: 
iBrosche. Odenkirchen fc GeffertI 'l999\ Moreover the im- 
pressive growth of computing facilities allow an accurate 
analysis of the uncertainties (a particularly critical point for 
this kind of application, see IHW for a discussion) by means 
of extensive Montecarlo simulations. 

In this paper I review the use of tidal radii of globu- 
lar clusters in the light of the more recent theoretical and 
observational result and I check if the mass estimates from 
this technique are consistent with "kinematical" estimates, 
dealing in particular with the mass profile at large Roc, e.g. 
the most interesting range (Sect. 2). The main conclusions 
and future prospects are summarized in Sect. 3. 



2 TIDAL RADII AS MASS PROBES 

A general formula for the tidal radius (rt) of a globular clus- 
ter of mass rwe orbiting around a gal axy of mass Mg{R) is 
jKindll962l:rOh. Lin fc Aarsethlll99l IHW, W99): 



rt 



k 



f{e)MG{Rp) 



Rr> 



(1) 



where Rp is the peri-galactic distance, /(e) is a function 
depending on the form of the galactic potential and on the 
eccentricity of the cluster orb it e and A; is a factor, firstly 
introduced bv lKeenanI l)l981allb|) to account for the elonga- 
tion of the limiting tidal surface along the line between the 
cluster center and the galactic center. According to this au- 
thor, IHW and W99 assume A; = |. From now on we drop, 
for brevity, the explicit indication of the dependence of mass 
on distance: in any case Mq shall be intended as the mass 
enclosed with a given galactocentric distance R (the corre- 
ponding index is also dropped, from now on R mea ns Roc)- 

The dependence on e was introduced bv lKind (^62) by 
reformulating the equation for the instantaneous tidal radius 
of Iron Hocrnor ( 1956) in terms of the perigalactic distance 
(instead of the instantaneous distance R). This choice was 
motivated by the fact that the typical orbital period (P) of 
globulars was expected to be much shorter than the internal 
relaxation time (tipically quantified by the half-mass relax- 
ation time trh)- In this case two-body relaxation is unable 
to keep the external structure of the cluster at pace with the 
changing galactic tidal field. In other words, the tidal force 
at perigalacticon truncates the cluster at the corresponding 
radius and the internal relaxation is too slow to restore a 
larger li miting radius before the next perigalactic passage 
(see also lOh. Lin fc Aarsethlll99,'Th . 

However recent comparisons between predicted and 
observed tidal radii that included all the information 
on clus ter orbits (as derived fro m meas ured proper 
motions iMeziane fc ColinI Il99d: J (^4epkir ch en et al.l Il99'^ 
iBrosche. Odenkirchen fc Gefferd Il999^ showed that ob- 
served tidal radii are larger than one would expect i f they 
were fixed at the perigalactic point. IMeziane fc ColinI l|l99(il l 
suggest that the actual tidal radius depends on the orbital 
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Figure 1. The logarithm of the ratio between the orbital period 
(P) and the half-mass relaxation time is plotted versus the galac- 
tocentric distance for all t he Galactic globulars for which orbita l 
p arameter are available in lDinescu. Girard fc van Altenal ^199g^ 
or lDinescu et all 1200 ll) . The horizontal line marks the threshold 
P/trh = li the vertical line marks Rqc = 35 Kpc. The names of 
the clusters with Rqc > 35 Kpc are also reported. 



phase, while IBrosche. Odenkirchen fc Gefferj lll999t) argue 
that a suitable average along the orbit defines a much more 
proper tidal radius with r espect to the perigalactic value. 

[Oh, Lin & Aarsptli lll992ll and Oh, Lin & Aarset^ 

11995f) studied in detail the properties of tidal radii of glob- 
ular clusters by means of N-body simulations. They found 
that while the tidal limit of clusters with P/trh << 1 is 
effectively set at perigalacticon, in clusters with P/trh ~ 1 
internal relaxation efi'ects are able to repopulate their exter- 
nal regions after the perigalacticon passage and their tidal 
radius may be comparable to the value expected at apogalac- 
ticon. 

In Fig. 1 the logarithm of the Qt = P/trh ra- 
tio is plotted versus the galactocentric distance for all 
the clusters for which orbital parameters from measured 
proper motions are available (from the computations b y 
iDinescu. Girard fc van Altenal Il999l : iDinescu et alJ 20W|V 
The re laxation times are from the 2003 version of the Harris! 
il996l) catalogue, which is the source of all the globular clus- 
ter data used in this study if not otherwise stated. There is a 
clear correlation between Qt and R suggesting that beyond 
7? ~ 35 kpc clusters are very likely to have Qt <'. 1, i.e., 
according to the results bv lOh. Lin fc AarsethI lll995l) . their 
ft does not reflect the tidal force at perigalacticon. 

On the basis of all the above considerations I make the 
assumption that the observed tidal radii of distant globu- 
lar clusters (e.g., those with R ^ 35 kpc) are probes of 
the Galactic tidal force at their present position. In other 
words I take the galactocentric distance Rp ^ R Ra 
as a reasonable approximation of the effective orbital ra- 
dius at which the tidal radius of clusters with Qt > 1 
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is fixed (IMeziane & Colin 


119961: iMevlan & Heesid Il997l: 


lOdenkirchen et alJ Il997l 




Brosche. Odenkirchen & GefFerd 


Il99fll: lOh. Lin fc AarsethI 


Il995h. The assumption limits the 



range that can be probed to the outer halo of the Galaxy. 
On the other hand this is the most interesting range since 
the contribution of barionic matter to the mass budget is 
expected to be negligible in this region. Furthermore, clus- 
ters orbiting in such remote regions are less likely to have 
their structure changed by close encounters with the galactic 
disk, bulge or with large molecular clouds, and are there- 
fore more reliable probes of the Galactic tidal field (see 
iMevlan'fc Hcgeic 1997, for an extensive review of the mech- 
anisms affecting the structure of GCs). 

Following the above assumption, the following form of 
the theoretical tidal radius is adopted: 



rt 



k 



2Mg 



R 



(2) 



e.g, Eq. 1 with /(e) = 2 and R instead of Rp. In the 
context of a logarithmic potential this is the exact fo rmula 
for the tidal radius of a cluster on a circular orbit iKind 

Imp)- 

Solving with respect to Ma- 



MG = ifc'me(^) . (3) 

W81, IHW and lOh. Lin fc AarsethI ^1995^ derived /(e) 
for eccentric orbits in a logarithmic potential. It is important 
to note that the adoption of /(e) — fiog{e), instead of the 
instantaneous value /(e) — fiog{& = 0) adopted here, has a 
modest impac t on the final Mg es t imate s. The adoption of 
the fiog{e) by lOh. Lin fc AarsethI il99 fJl changes the final 
Mg estimate by a factor ^ 2 for e 0.825 and by a factor 
^ 4 for e ^ 0.935 with respect to the /(e) = 2 case adopted 
here. In any case the derived Ma is larger than what ob- 
tained with Eq. 3. In analogy with the analysis that will be 
described in Sect. 2.2 I performed an extensive set of Mon- 
tecarlo experiments using the fiog{e) version of Eq. 1 and 
exploring the whole range of possible orbital eccentricities. 
The resulting mass (and mass profile) estimates are fully 
consistent with the results obtained from Eq. 3. Hence, the 
approximations involved in the adopted appr oach appear 
fully adequate for the present purpose (see also llbata et alJ 

2.1 The observational side 

Tidal radii of GCs are not, in general, directely observ- 
able quantities. Fitting King models (King 1962) to surface 
brightess profile of globular clusters one obtains an estimate 
of the core radius (rc) and of the concentration parameter 
C = log (■^); rt is derived from these fitted parameters. 
Hence, the best estimates of rt may be obtained from bright 
clusters (e.g., providing high signal-to- noise data for the fit) 
whose surface brightness profile is reliably measured over the 
largest possible radial range (thus limiting at a minimum the 
extrapolation to the actual tidal limit). 

Unfortunately the range Roc > 35 kpc is mainly pop- 
ulated by sparse, low luminosity clusters whose tidal radii 
and integrated magnitudes are quite uncertain. The suitable 
galactic clusters in the relevant range of galactocentric dis- 
tances are: Pal 15, NGC 7006, Pyxis, Pal 14, NGC 2419, 
Eridanus, Pal 3, Pal 4 and AM-1. Pal 15 and Pyxis have 



been excluded since they are affected by significant amount 
of interstellar extinction {Ay > 0.6 mag, that implies a 
larger uncertainty in the estimate of rric). Note however that 
the inclusion of these clusters does not change in any way 
the final results presented below. Of the remaining clusters 
the best suited for the present analysis are NGC 7006 and 
NGC 2419, two bright and well studied clusters. In particu- 
lar NGC 2419 is the 4-th most luminous cluster of the whole 
Galaxy {My ~ —9.58) and its surface brightness profile have 
been reliably measured out to the 85 % of the deduced tidal 
radius, i.e. the level of extrapolation is quite modest. Finally 
NGC 2419 is the only cluster of the considered set for which 
a direct estimate of the mass-to-light ratio {M/L, a funda- 
ment al ingredient to derive rric) is available {M/Lv = 1.2, 
from iPrvor fc Mevlanlll993l) . 

In conclusion, NGC 7006 and NGC 2419 provide by far 
the most robust and less uncertain mass probes, the remain- 
ing clusters are retained just for consistency check. 

2.1.1 Montecarlo simulations 

Mg is estimated from each cluster using Eq. 3. To deal 
with uncertainties I obtained 10000 independent Mg esti- 
mates for each cluster by extracting at random (from suit- 
able distributions described below) the following quantities: 
the observed V-band distance modulus jj,v, the apparent in- 
tegrated V magnitude Vt, the tidal radius in arcmin rt, the 
V-band mass-to- light ratio in solar units M/L and Keenan's 
k factor (see below). The color excess [E{B — V)] is kept 
fixed, since it is small (less than ~ 0.1 mag) in all of the 
considered cases. For each set of extracted parameter the 
following items are computed: the distance from the Sun 
and the Galactic Center from fiv, E{B — V) and the galac- 
tic coordinates, the linear tidal radius from its angular value 
and distance, the mass of the cluster from Vt, jJ.v and M/L, 
and finally, from R and rt in kpc and rric in solar masses, 
the Galactic mass enclosed within R. All the distributions 
are chosen to (conservatively) cover the whole range that is 
compatible with the adopted uncertainties of each param- 
eter, hence the final 10000 Mg estimates cover the whole 
range allowed by taking into account all the possible sources 
of error. Finally the median of the 10000 estimates is com- 
puted as well as the range in Mg including the 90 % of the 
derived Mg estimates. 



2.1.2 Uncertainties on input parameters. 

To have a closer look to the details of the simulations, the 
assumptions on input parameters are shortly described be- 
low. 

• k factor. Keenan's factor (fc = |) provides an av- 
erage correction for the non-spherical shape of the limit- 
ing tidal surface. To account for the possible cluster-to- 
cluster variation of this parameter, fc is extracted from a 
uniform distribution in th e range 0.5 ^ fc ^ 1.0 (see also 
iHeggie fc Ramamanill995h . 

• fly. To account for both the measurement errors and 
the uncertai nties still affec ting the distance scale of globular 
clusters (see lCacciarll999l . and references therein) values 
were extracted from a gaussian dist ribution with mean equal 
to the listed by iHarrij il996ll and standard deviation 
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= 0.1 mag for NGC 2419 and NGC 7006 and cr^ = 0.2 
mag for the remaining, less extensively studied, clusters. 

• Vt values were extracted from a gaussian distribution 
with mean equal to the Vt listed bv THarris ( 199£) and stan- 
dard deviation = 0.1 mag for NGC 2419 and NGC 7006 
and av ~ 0.2 mag for the remaining clusters. Note that 
the measure of this parameter is particularly critical for 
low-density and low-brightness clusters, like the majority 
of those considered here. 

• M/L. The observed Funda mental Plane of GCs 
JPiorgovsy I1995I: iBellazzini 1 19981) implies that the M/L 
ratio of globular clusters is constant to within a factor 
of a few, a pproximately compatib le with the measure- 
ment errors. IPrvor fc Mevlai] iim^t find 0.5 ? M/L ~ 4 
in goo d agreement with lMaiidiishev. Soassova fc Staneval 
il99ll) . Here the M/L values are extracted from a gaus- 
sian distribution with mean M/L = 1.2 (e.g., the esti- 
m ated M/L of NGC 2419 and the aver age M/L derived 
bv lMandushev. Soassova fc StanevalliggiT) and ctm/l = 0.4, 
with the further constraint M/L > 0.5. This assumption 
ensures that the range of observed M/L of GCs is fully ex- 
plored. 

• rt- I searched the literature to find estimates of rt for 
th e considered clusters that ma y supersede those reported 
bvlHarrij (Il996l) , all drawn from lTraeer. King fc Diorgovskil 
^1995^ . The only (partially) succ essful case was Pal 14, 
for wh ich I adopt the estimate bv lHarris fc van den Berghl 
il984l) whose density p rofile is marginally more extende d 
with respect to that bv lTrager. King fc Diorgovskil il995f) . 
The adopted rt values are extracted from a gaussian distri- 
bution with the mean equal to the listed values and with 
standard deviation in the range ar — O.lr* — 0.3rt, depend- 
ing on the quality and extension of the available surface 
brightness profile. In particular Or ~ 0-lrt for NGC 2419, 
ar = 0.2rt for NGC 7006 and ar = O.Srt for the remaining 
clusters. 

The Referee correctly pointed out that the estimates 
of tidal radii may be also plagued by systematics. For ex- 
ample the adoption of different models (e.g. IWilson 19751) 
to fit the surface brightness profiles of globulars may lead 
to obtain significantly larger limiti ng radii than what es- 
timat ed with King's models (see iMcLaughlin fc MevlanI 
I2OO3! . for an application). Moreover, some theoretical stud- 
ies also suggest that realistic models of globular clusters 
in the Galactic tidal field may be sligh tly more spatially 
extended th an King' s models (see, e.g. iKashlinskvl Il988t 
[Hcggic & Ramaman|[l995). To explore the effect of sys- 
tematics that may change the observationally estimated 
Vt values up to a factor ~ 2 I repeated the analysis de- 
scribed above for NGC 2419 ( rt = 8.74', according to 
iTrager. King fc Diorgo^ Il995fl assuming ar = 0.5rt in- 
stead of ar = O.lrt. To preserve compatibility with the ob- 
served profile, which reaches r ~ 7.5', I forced the randomly 
extracted tidal radii to the range rt > 6'. Hence the final rt 
of the Montecarlo simulation are in the range 6' < rt < 20'. 
Note that at r J5 6' the observed profile is rapidly falling 
and the surface brightness is Ev > 28 mag/arcsec^, thus 
it is quite unlikely that the actual limiting radius is much 
larger than 10'. With these new assumptions the typical 
uncertainty on the final Mg estimates grows from ~ 75% 
(for the case ar = O.lrt) to ~ 110% while the derived me- 



dian Mg is practically unchanged. Therefore the inclusion 
of these possible systematic errors in the uncertainty budged 
doesn't seriously affect the main conclusions of this paper 
(see Sect. 2.2), at least for what concern the cluster that pro- 
vide the most interesting constraint on Mg, e.g. NGC 2419. 

2.2 The mass of the Galaxy within 35 < Rgc < 100 
kpc. 

Fig. 2 reports the median Mg and R of the 10000 simulation 
carried on for each considered cluster (filled circles). The 
error bars enclose 90 % of the derived estimates (±45 %). 
The most reUable points (NGC 7006 and NGC 2419) are 
indicated by larger symbols. 

Other mass estimates, more recent than those consid- 
ered by Z99, are also reported (sec legend). The only pre- 
99 estimate reported is that by iKochanek (1996), who pro- 
vides the more extensive treatment of the problem at that 
epoch. The filled triangle with an arrow, labeled L-B99 in 
the legend, provides a sensible upper limit to the total mass 
of the Galaxy since it is the estimate of th e mass of the 
Local Group obtained by [Lvnd en-Bell 1^9^ using the Lo- 
cal Group timing technique. The reported error bars have 
etherogeneous meanings (1 — ct errors, 90 % c.L, etc., as 
provided by the authors) and hence they are not directly 
comparable. The mass profiles of isothermal spheres with 
Vc = 180, 220, 260 km/s are also plotted for reference. They 
are intended to allow the comparison of the estimates at 
large R with constraints provided by the rotational veloc- 
ity of the HI disk at ii < 20 kpc, taking into account 
the wh ole range of possible uncertainty on the Galactic 
Vc (s ee iFitch fc Tremainelll99ll : iPedreros. Anguita fc Mazal 
I2OO2L Z99, and references therein). 

From the inspection of Fig 2 the following main conclu- 
sions can be drawn: 

(i) All the reported estimates based on kinematical data 
are in agreement, within the uncertainties, with the reported 
mass profiles and with the results summarized by Z99. 

(ii) All the estimates based on the tidal radii of globular 
clusters are consistent with the reported mass profile, within 
the (large) uncertainties. 

(iii) The clusters that provide the most reliable constraint 
with the adopted technique (e.g., NGC 7006 and NGC 2419) 
give estimates of the enclosed Galactic mass in excellent 
agreement with the reported profiles and with all other 
estimates. In particular, from NGC 7006 I obtain Mg (39 
kpc) = 0.381J J X 10^^ Mq, and from NGC 2419, Mg (92 
kpc) = l.Stli'x 10^2 Mq. 

Hence, the present application of the tidal radii tech- 
nique provide an independent validation of the standard 
framework for the mass of the Milky Way as it emerges from 
the analysis by Z99 and from more recent studies. There is 
a quite remarkable general agreement among all the consid- 
ered estimates, indicating that (a) the mass of the Galaxy 
appears to grow with galactocentric distance at least up to 
R > 100 Kpc and (b) the total mass of the Galaxy is larger 
than ~ 2 X 10^^ Af©- This implies that the V-band mass to 
light ratio of the Milky Way spans a range 20 < M/L < 100 
in the radial range 30 kpc < i? < 200 kpc (see Fig. 3 by 
Z99). Hence, in the classical Newtonian theory, all the avail- 
able observational constraints consistently point to the con- 
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Figure 2. The Mass profile of the Milky Way halo (for Rqc > 35 kpc) from various recent estimates based on kinematical data (see legend 
in the lower left corner) and from the tidal radii of the selected globular clusters (filled circles). Large filled circles are the best observed 
globulars, providing (by far) the stronger constraint on Mq with respect to the other clusters. The error bars for globulars enclose the 
90 % c.l. range. The mass profiles of isothermal sphe res with Vc = 180, 220, 260 km/s are also rep orted for ref erence. Acronyms: WE99 = 
IWilkinson fc Evan j il99fl^: K96 =lKqchanekl Jl99(J) : S CB03 = fSakamoto. Chiba Beerj ^200^^) : PAM02 = IPedreros. Anguita Mazal 
^2002^ : vdM03 = lvan der Marel et alj 12003) : L-B99 = lLvnden-Bcll 1 1999) . mass of the Local Group. The arrow indicates that the last 
one is an upper limit for the total mass of the Galaxy. 



elusion that the Milky Way is surrounded by a huge Dark 
Matter halo. 



3 CONCLUSIONS 

The tidal radii of galactic globular clusters have been 
used to obtain estimates of the enclosed mass of the 
Galaxy in the range of galactocentric distances 35 kpc 
^ Rgc ^ 100 kpc, under the assumption that, for 



such remote clusters, Rgc is a reasonable approxima- 
tion of t he orbital radius at which their tidal limit is 
imposed (Meziane & Colin '199^; Odenkirchen et al. ''1997|; 
l^oscho, Odenkirchen & Gcffcrt 1999; Oh, Lin & Aarss^ 
Il995fl . The adopted technique provides an estimate of the 
enclosed galactic mass that is independent of kinematical 
data that, on the other hand, provides the basis of all other 
existing estimates (see Sect. 1). Therefore, while the asso- 
ciated uncertainties are still quite large, the present appli- 
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cation provide (at least) an interesting consistency check of 
the existing kinematical estimates of Mq- 

The estimates obtained from tidal radii are fully con- 
sistent with results from other authors and other methods. 
All the available constraints are consistent with the mass 
profile of an isothermal sphere with Vc = 220 ± 40 km/s. 
Hence, the present analysis provides independent support 
to the fact that the mass of the Galaxy grows with R out to 
large distances from the Galactic Center and that the mass 
enclosed within J? ~ 90 kpc is Mg 10^^ Mq . 

The present analysis indicates that tidal radii of remote 
globular clusters may be better probes of the galactic poten- 
tial than previously believed (IlfW) . A detailed and compre- 
hensive theoretical analysis it is now at hand with rea listic 
N-body simulations (e.g., with GRAPE. lMakindll99^ and 
can provide much sounder and solid basis to the technique. 
At the same time, wide-field cameras mounted on large tele- 
scopes may provide the opportunity to obtain more extended 
density profiles, based on star counts, hence reducing the ob- 
servational uncertainty. These advancements may ultimately 
lead to a fully reliable additional technique to probe the mass 
profile of our Galaxy, nicely independent of and complemen- 
tary to the usual kinematical methods. 
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